A New Limit on the Flux of Cosmic Antihelium* 
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A very sensitive search for cosmic-ray antihelium was performed using data obtained from three 
scientific flights of BESS magnetic rigidity spectrometer. We have not observed any antihelium; this 
places a model-independent upper limit (95 % C.L.) on the antihelium flux of 6 x 10~ 4 m -2 sr _1 s _1 
at the top of the atmosphere in the rigidity region 1 to 16 GV, after correcting for the estimated 
interaction loss of antihelium in the air and in the instrument. The corresponding upper limit on the 
He/He flux ratio is 3.1 xlO -6 , 30 times more stringent than the limits obtained in similar rigidity 
regions with magnetic spectrometers previous to BESS. 

PACS numbers: 98.80.Cq, 98.70.Sa, 98.80.Bp, 98.90.+S 



Cosmic-ray observations provide most direct evidence 
for our Galaxy being composed mostly by baryons. This 
baryon-antibaryon asymmetry can be global in the Uni- 
verse, being created in the very early Universe due to the 
violations of CP and of baryon-number. However, de- 
pending on the nature of CP violation, baryon-symmetric 
models are conceivable in which the Universe is sep- 
arated into an equal number of matter- and antimatter- 
domains. Whereas 7-ray observations place strong limi- 
tations on the antimatter in our Galaxy and in the local 
cluster of galaxies, the domain structure could still ex- 
ist beyond this scale. Although antihelium might be in 
principle produced in cosmic-ray interactions, their con- 
tribution to the He/He flux ratio is expected to be much 
smaller than 1CP 12 |2j. Detecting antihelium at a level 
higher than this could therefore provide the evidence of 
antimatter domains or of other exotic phenomena such 
as superconducting strings in our Galaxy Q . For further 
discussion of astrophysical considerations regarding the 
search for antihelium see the references [f|,|]. We report 
in this letter (for more detail see ||) a sensitive search 
for antihelium using the data from the '93, '94, and '95 
flights of BESS detector, and provide model-independent 
upper limits on the absolute antihelium flux as well as on 
the He/He flux ratio. A limit on the flux ratio from an 
early analysis of '95 data is published elsewhere 

Fig. [I] shows front- and side-views of the BESS '95 in- 
strument. The cylindrical configuration provides a wide 
tracking region and an acceptance of up to 0.32 m 2 sr 
depending on the off-line fiducial cuts. From inside to 
outside, it includes a jet-type drift (JET) chamber, in- 
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ner drift chambers (IDCs), a superconducting solenoid, 
outer drift chambers (ODCs), and a time of flight (TOF) 
hodoscope. The solenoid produces a magnetic field of 1 
T with an uniformity of ±15% inside the bore. 

The JET chamber 0, as the key tracking detector, 
measures up to 24 points per track thrcc-dimcnsionally, 
each with the resolution |gj of 200 fim in r<p plane and of 2 
cm in z-position (a cylindrical coordinate (rcj)z) is deter- 
mined by defining the magnetic field direction as z-axis) . 
Each of the IDCs and ODCs consists of two 12-mm- 
thick drift layers which are divided into 50- mm- wide cells, 
and measures r<j> position each with 200 /xm resolution. 
These redundant and continuous position measurements 
with ODCs, IDCs, and JET chamber, all equipped with 
multi-hit capacity, make it possible to recognize multi- 
track events and tracks having interactions and scatter- 
ings, thus minimizing the background originating from 
the interactions. The r</>-tracking in the central region 
is performed by fitting hit points in the JET chamber 
as well as in the IDCs. This results in a maximum de- 
tectable rigidity (R) of 200 GV and a typical rigidity 
resolution of 0.5 % at 1 GV. Using vernier pads with a 
cycle of 10 cm for IDCs and 12 cm for ODCs, both cham- 
bers can measure the z position, modulo the cycle, with 
300 //m resolution. Combining the IDC and JET infor- 
mation, z-position of the track can be determined with a 
precision of 300 /im. 

The TOF hodoscope, placed at a radius of 65 cm, con- 
sists of eight upper and twelve lower plastic scintillators 
(Bicron 404), each of which has a dimension of 95 cm 
x 10 cm x 2 cm. The light signals from the scintilla- 
tor are guided through acrylic light-guides and reach the 
photomultiplier-tubes (PMTs) attached on both ends. 
The timing and amplitude of the PMT signals are mea- 
sured to determine the time of flight and dE/dx of inci- 
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dent particles with resolution of 110 ps and 10 %, respec- 
tively. For '93 and '94 instruments, the TOF hodoscope 
consisted of four upper and six lower plastic scintillators, 
each of which had a dimension of 110 cm x 20 cm x 2 
cm. The timing and dE/dx resolutions was 280 ps and 15 
%, respectively. A detailed description of '93 instrument 
has been published elsewhere 0|. 

These detector-components along with electronics are 
enclosed by a 2.5-mm-thick aluminum vessel that keeps 
the inside pressure at 1020 mb during the flights. The 
total material-thickness along the particle path is 17.8 
g/cm 2 , of which 8.6 in front of the central tracker, 0.4 
between IDC and JET, and 0.14 in the JET chamber 
gas, all in g/cm 2 . 

The trigger system is designed to detect negatively 
charged particles (p, antihclium) efficiently while sam- 
pling the proton and helium events in an unbiased man- 
ner. This system is composed of the "TO trigger" and 
the "track trigger". The TO trigger initiates the data 
gathering based upon a simple coincidence of the top 
and bottom scintillators with thresholds set at 1 /3 of the 
pulse height from a minimum ionizing particle. The track 
trigger consists of two sequential processes. The first is 
a pattern-selection process to reject the null-track and 
multi-track events by using the hit-pattern of the IDC 
and ODC cells. If an event passes the pattern-selection, 
the second process performs a rough rigidity-selection by 
using the cell hit-pattern. Irrespective of the track trigger 
condition, a fraction of TO triggers is recorded to provide 
an "unbiased trigger sample" , which is used to deter- 
mine the efficiency of the track trigger selections. The 
data were collected in balloon-flights performed on July 
26 '93, July 31 '94, and July 25 '95, all launched from 
Lynn Lake, Manitoba, Canada. During each level flight, 
the floating altitude was 36.5 km (residual atmospheric 
depth of 5 g/cm 2 ). The total live time of the three flights 
was 84446 sec. 

The following off-line selections are applied for both 



negative and positive curvature events, (i) One and only 
one track should be found in the JET chamber, (ii) The 
track should be fully contained in the central six-column 
of the JET chamber (see fig[l]), where the resolution of the 
rigidity-measurement is optimum. This fiducial cut de- 
termines the geometrical acceptance for the present anal- 
ysis of 0.35 m 2 sr for '93 and '94, and 0.28 m 2 sr for '95. 
The following cuts are then applied to ensure the quality 
of the track and the timing measurement. (1) The fit- 
ted r</>-track should contain at least 14 hits in the JET 
chamber and at least one hit in each of the upper two 
and lower two IDC layers. (2) The reduced chi-square of 
the r</>-track-fitting has to be less than 3. (3) The extrap- 
olated track should cross the fiducial region of the TOF 
scintillators (\z\ < 49.5 cm for '93 and '94, \z\ < 47.0 cm 
for '95). (5) The z-position (ztof) determined by the 
left-right time difference measured by the PMTs should 
match the z-impact point of the extrapolated track at 
the TOF counter within 10 cm (|ztof — z\ < 10 cm). 

The particle velocity, i.e., /?(= v/c) is determined from 
the time of flight and the path-length. Upward-moving 
particles are completely separated from downward- 
moving ones by the sign of (3 due to the excellent time res- 
olution. We reject all upward-moving (albedo) particles 
at this stage, and limit further analysis to the downward- 
moving particles. Fig. ^ (a) and (b) show plots of dE/dx 
versus \R\ for the top and bottom scintillators, respec- 
tively, for a sample of '95 data after applying all of the 
above mentioned cuts. In order to select helium and anti- 
hclium, we require that events must reside in the ll dE/dx 
bands" shown in Fig. || (a) and (b). These dE/dx cuts 
reject most of singly-charged particle events (p, p, e ), 
while preserving a high efficiency (95 %) for helium. Fig. 
H shows plot of 1/(3 versus \R\ for a sample of '95 data 
after these dE/dx cuts, where clear lines of 4 He and 3 He 
are visible. We further require that antihclium as well as 
helium have f3 and rigidity values inside the band shown 
in Fig. [| The efficiency of the /3-band cut is 99 % for 
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FIG. 1. Cross-sectional front- and side-views of the BESS '95 instrument. 
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FIG. 2. Scatter plot of AE/dx vs. absolute rigidity for (a) 
the top and (b) bottom layers of the TOF hodoscope ('95). 
The solid lines show the d_E/da; band cut positions. 

helium. The overall efficiency of the off-line helium se- 
lection (from cut (1) to the /3-band cut), is 63 %, 70 %, 
and 78 %, for '93, '94, and '95, respectively, and 71 % in 
average (weighted by number of events). 

All selections up to this level use \R\ instead of R, 
i.e. do not discriminate the positive and negative rigidity 
(helium and antihelium). Fig. || shows the 1/R distribu- 
tion of the events which survived all the cuts mentioned 
above. Those events in the negative region are appar- 
ently the spillover from the positive region due to the 
finite 1 /R resolution. There exist no events to the left of 
-0.0625 GV" 1 , which corresponds to the rigidity of -16 
GV. To define the lowest end of the rigidity region, we 
note that the number of helium in Fig. |] drops sharply 
at 1.4 GV -1 , which corresponds to 1 GV at the top of 
the atmosphere (TOA) after correcting for the energy 
losses in the instrument and in the air. Therefore, we 
conclude that we observed no antihelium in the rigidity 
region from 1 to 16 GV at TOA. 

The corresponding limit on the He/He flux ratio at 
TOA can be obtained by dividing the upper limit on 
the number of antihelium (Nn m ) by the total number 
of helium (N), which can be expressed as the integral; 
N = X^i[ n i/( e pat£sei??)i]/f- The summation is over the 60 
bins which divide the rigidity region (1 to 16 GV), nj is 
is the number of helium detected in the unbiased trigger 
sample falling in the i-th rigidity bin, e pa t and e sc i are the 
efficiencies, for the helium, of the track-pattern- and the 
off-line selections, respectively, and rj is the probability of 
the helium surviving through the air and the instrument, 
all at the i-th rigidity bin. The pre-set sampling factor f 
of the unbiased trigger was 1/40 for '93, 1/15 for the first 
half and 1/30 for the last half of '94, and 1/20 for '95. 
The efficiencies e pat as well as e so i can be directly deter- 
mined by using the unbiased trigger sample. To calculate 
the surviving probability 77, we utilize the "inelastic" fic[ | 
cross sections of helium incident on various target nuclei, 
i.e. <r(He, A t ). Since the data on <r(He, A t ) at GeV ener- 
gies is relatively sparse, we obtain cr(He, At) by starting 
with a(p, A t ), given as a function of incident energy [jn], 
and by scaling it using "the hard sphere with overlap" 
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FIG. 3. Scatter plot of 1//3 vs. absolute rigidity. The solid 
lines show the cut positions. ('95) 
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FIG. 4. 1/rigidity distribution of the events which survive 
all selections. 



model of nuclear interaction |12|, i.e. 

a(A u A t ) ex (Al /3 + Al /3 - 0.71 x {A7 1/3 + A; 1/3 )) 2 (1) 

where A\ and At are the atomic weight of incident and 
target nuclei. This model is known to reproduce data 
on nuclear interactions for various combinations of A\ 
and At including the proton |jl3). The resultant c(He, 
At) agrees (within an accuracy of 10 %) with measured 
helium cross sections on targets such as C or Al. This 
also justifies our use of the same formula ([j]) to estimate 
the unmeasured c(He, A t ) cross sections from the a(p, 
A t ) data @. 

When summed over the rigidity bins, the number of 
detected helium events divided by the sampling factors 
(E«i/f) are 2 -3 x 10 5 , 8.7 x 10 5 , and 8.6 x 10 5 for '93, 
'94, and '95, respectively, and the total number of helium 
N is 5.1 x 10 6 . _ 

The upper limit on the number of antihelium (N) de- 
pends on its unknown rigidity spectrum, since the ef- 
ficiencies are functions of rigidity. In order to obtain 
the most conservative limit, irrespective of the antihe- 
lium spectrum, we define N Um = 3.0/[e pat e r ige se r?7]min, 
where 3.0 is the 95 % confidence level limit for zero- 
observed events, and [e pa te r i g e S ci?7]min is the minimum 
value of the antihelium efficiency-product in the rigidity 
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FIG. 5. The resultant upper limit of He/He flux ratio of 
the present paper together with previous limits. 

region. In calculating [e pa te r ige se i7?] m in, we take e se i = e sc i 
and e pat = e pa t for each rigidity bin, since the pattern- 
selection in the track trigger as well as the off-line se- 
lection does not depend on the sign of the track rigidity. 
The efficiency of the rigidity-selection in the track trigger 
e r i g can be directly obtained by using events with a neg- 
ative rigidity track in the unbiased trigger sample, and 
is a monotonic function of the rigidity, having values of 
0.95 at 1 GV and 0.85 at 16 GV. The efficiency prod- 
uct (e pa te r i g e S ei77) thus does not vary significantly (0.19 
at 0.7 GV, 0.23 at 2 GV, and 0.20 at 16 GV). We take 
the minimum value of 0.19 to calculate Ni; m . 

The resultant 95% C.L. upper limit on the He/He flux 
ratio at TOA is 3.1 x 10~ 6 in the rigidity region from 1 to 
16 GV. The upper limit on the antihclium flux integrated 
over the rigidity region is 6 x 10~ 4 m _2 sr -1 8 -1 . It should 
be noted that these upper limits are very conservative 
and are valid for any hypothetical antihelium spectrum. 
Our limit on the flux ratio obtained by this work is shown 
in Fig. U together with the previous ones |^,^5). Our 
flux ratio is a factor 7 improvement over the limit of 
Buffington et al. who looked for the annihilation signal 
of low rigidity (1 ~ 2 GV) antihclium in a spark chamber 
calorimeter, and is 30 times more stringent than the limit 
of Golden et al. who covered a rigidity region similar 
to ours by using a previous generation magnetic rigidity 
spectrometer. The large acceptance of BESS, which is 
designed ||l6|| specially for the antimatter and antiproton 
searches, made it possible to set this limit. 
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